%to plot CIS locations
[map,maplegend] = gtopo30('c:\manoj\download\topo\E060N40.DEM',1,[20,22],[75,78]);
[a,b] = xlsread('c:\manoj\temp\cis-loc.xls');
sstn(1:length(a),1) = ' ';
stn = char(b(:,1));
stn = [sstn stn];
he = ltln2val(map,maplegend,a(:,1),a(:,2));
plot3m(a(:,1),a(:,2),he+100,'k.','MarkerSize',20);
textm(a(:,1),a(:,2),he+100,stn);